Revealing modified gravity signal in matter and halo hierarchical clustering

We use a set of N-body simulations employing a modified gravity (MG) model with Vainshtein screening to study matter and halo hierarchical clustering. As test-case scenarios we consider two normal branch Dvali-Gabadadze-Porrati (nDGP) gravity models with mild and strong growth rate enhancement. We study higher-order correlation functions $\xi_n(R)$ up to $n=9$ and associated hierarchical amplitudes $S_n(R)\equiv\xi_n(R)/\sigma(R)^{2n-2}$. We find that the matter PDFs are strongly affected by the fifth-force on scales up to $50h^{-1}$Mpc, and the deviations from GR are maximised at $z=0$. For reduced cumulants $S_n$, we find that at small scales $R\leq10h^{-1}$Mpc the MG is characterised by lower values, with the deviation growing from $7\%$ in the reduced skewness up to even $40\%$ in $S_5$. To study the halo clustering we use a simple abundance matching and divide haloes into thee fixed number density samples. The halo two-point functions are weakly affected, with a relative boost of the order of a few percent appearing only at the smallest pair separations ($r\leq 5h^{-1}$Mpc). In contrast, we find a strong MG signal in $S_n(R)$'s, which are enhanced compared to GR. The strong model exhibits a $>3\sigma$ level signal at various scales for all halo samples and in all cumulants. In this context, we find that the reduced kurtosis to be an especially promising cosmological probe of MG. Even the mild nDGP model leaves a $3\sigma$ imprint at small scales $R\leq3h^{-1}$Mpc, while the stronger model deviates from a GR-signature at nearly all scales with a significance of $>5\sigma$. Since the signal is persistent in all halo samples and over a range of scales, we advocate that the reduced kurtosis estimated from galaxy catalogues can potentially constitute a strong MG-model discriminatory as well as GR self-consistency test.Comment: 19 pages, 11 figures, comments are welcom

Improved analytical modeling of the non-linear power spectrum in modified gravity cosmologies

Reliable analytical modeling of the non-linear power spectrum (PS) of matter perturbations is among the chief pre-requisites for cosmological analyses from the largest sky surveys. This is especially true for the models that extend the standard general-relativity paradigm by adding the fifth force, where numerical simulations can be prohibitively expensive. Here we present a method for building accurate PS models for two modified gravity (MG) variants: namely the Hu-Sawicki $f(R)$, and the normal branch of the Dvali-Gabadadze-Porrati (nDGP) braneworld. We start by modifying the standard halo model (HM) with respect to the baseline Lambda-Cold-Dark-Matter ($\Lambda$CDM) scenario, by using the HM components with specific MG extensions. We find that our $P(k)_{\text{HM}}$ retains 5% accuracy only up to mildly non-linear scales ($k \lesssim 0.3$ h/\,\mbox{Mpc}) when compared to PS from numerical simulations. At the same time, our HM prescription much more accurately captures the ratio $\Upsilon(k) = P(k)_{\text{MG}}/P(k)_{\Lambda \text{CDM}}$ up to non-linear scales. We show that using HM-derived $\Upsilon(k)$ together with a viable non-linear $\Lambda$CDM $P(k)$ prescription (such as HALOFIT), we render a much better and more accurate PS predictions in MG. The new approach yields considerably improved performance, with modeled $P(k)_{\text{MG}}$ being now accurate to within 5% all the way to non-linear scales of $k \lesssim 2.5-3$ h/\,\mbox{Mpc}. The magnitude of deviations from GR as fostered by these MG models is typically $\mathcal{O}(10\%)$ in these regimes. Therefore reaching 5% PS modeling is enough for forecasting constraints on modern-era cosmological observables

Uneven flows: On cosmic bulk flows, local observers, and gravity

Using N-body simulations we study the impact of various systematic effects on the bulk flow (BF) and the Cosmic Mach Number (CMN). We consider two types of systematics: those related to survey properties and those induced by observer's location in the Universe. In the former category we model sparse sampling, velocity errors, and survey incompleteness. In the latter, we consider Local Group (LG) analogue observers, placed in a specific location within the Cosmic Web, satisfying various observational criteria. We differentiate such LG observers from Copernican ones, who are at random locations. We report strong systematic effects on the measured BF and CMN induced by sparse sampling, velocity errors and radial incompleteness. For BF most of these effects exceed 10\% for scales $R\leq100h^{-1}$Mpc. For CMN some of these systematics can be catastrophically large ($>50\%$) also on bigger scales. Moreover, we find that the position of the observer in the Cosmic Web significantly affects the locally measured BF (CMN), with effects as large as $\sim20\%$ ($30\%)$ at $R\leq50h^{-1}$Mpc for a LG-like observer as compared to a random one. This effect is comparable to the sample variance. To highlight the importance of these systematics, we additionally study a model of modified gravity (MG) with $\sim15\%$ enhanced growth rate. We found that the systematic effects can mimic the modified gravity signal. The worst-case scenario is realized for a case of a LG-like observer, when the effects induced by local structures are degenerate with the enhanced growth rate fostered by MG. Our results indicate that dedicated constrained simulations and realistic mock galaxy catalogs will be absolutely necessary to fully benefit from the statistical power of the forthcoming peculiar velocity data from surveys such as TAIPAN, WALLABY, Cosmic Flows-4 and SKA.Comment: 20 pages, 9+2 figures, comments are welcome

A dark matter solution to the H0H_{0} and σ8\sigma_{8} tensions, and the integrated Sachs-Wolfe void anomaly

We consider a phenomenological model of dark matter with an equation-of-state that is negative and changing at late times. We show this couples the $H_{0}$ and $\sigma_{8}$ tensions, providing an explanation for both simultaneously, while also providing an explanation for the anomalously large integrated Sachs-Wolfe (ISW) effect from cosmic voids. Observations of high ISW from cosmic voids may therefore be evidence that dark matter plays a significant role in the $H_{0}$ and $\sigma_{8}$ tensions. We predict the ISW from cosmic voids to be a factor of ~ 2 greater in this model than what is expected from the standard model $\Lambda$CDM.Comment: 5 pages (+references) and 3 figures. Comments welcom

Caught in the cosmic web:Environmental effect on halo concentrations, shape, and spin

Using a set of high-resolution simulations we study the statistical correlation of dark matter halo properties with the large-scale environment. We consider halo populations split into four Cosmic Web (CW) elements: voids, walls, filaments, and nodes. For the first time we present a study of CW effects for halos covering six decades in mass: $10^{8}-10^{14}{h^{-1}{\rm M}_{\odot}}$. We find that the fraction of halos living in various web components is a strong function of mass, with the majority of $M>10^{12}{h^{-1}{\rm M}_{\odot}}$ halos living in filaments and nodes. Low mass halos are more equitably distributed in filaments, walls, and voids. For halo density profiles and formation times we find a universal mass threshold of $M_{th}\sim6\times10^{10}{h^{-1}{\rm M}_{\odot}}$ below which these properties vary with environment. Here, filament halos have the steepest concentration-mass relation, walls are close to the overall mean, and void halos have the flattest relation. This amounts to $c_{200}$ for filament and void halos that are respectively $14\%$ higher and $7\%$ lower than the mean at $M=2\times10^8{h^{-1}{\rm M}_{\odot}}$, with low-mass node halos being most likely splashed-back. We find double power-law fits that very well describe $c(M)$ for the four environments in the whole probed mass range. A complementary picture is found for the average formation times, with the mass-formation time relations following trends shown for the concentrations: the nodes halos being the oldest and void halo the youngest. The CW environmental effect is much weaker when studying the halo spin and shapes. The trends with halo mass is reversed: the small halos with $M<10^{10}{h^{-1}{\rm M}_{\odot}}$ seem to be unaffected by the CW environment. Some weak trends are visible for more massive void and walls halos, which, on average, are characterized by lower spin and higher triaxiality parameters.Comment: 18 pages, 9 figures, match the published version in Physical Review D eid. 06351